home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / zlaqgb.z / zlaqgb
Text File  |  1998-10-30  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZLLLLAAAAQQQQGGGGBBBB((((3333FFFF))))                                                          ZZZZLLLLAAAAQQQQGGGGBBBB((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZLAQGB - equilibrate a general M by N band matrix A with KL subdiagonals
  10.      and KU superdiagonals using the row and scaling factors in the vectors R
  11.      and C
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE ZLAQGB( M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND, AMAX,
  15.                         EQUED )
  16.  
  17.          CHARACTER      EQUED
  18.  
  19.          INTEGER        KL, KU, LDAB, M, N
  20.  
  21.          DOUBLE         PRECISION AMAX, COLCND, ROWCND
  22.  
  23.          DOUBLE         PRECISION C( * ), R( * )
  24.  
  25.          COMPLEX*16     AB( LDAB, * )
  26.  
  27. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  28.      ZLAQGB equilibrates a general M by N band matrix A with KL subdiagonals
  29.      and KU superdiagonals using the row and scaling factors in the vectors R
  30.      and C.
  31.  
  32.  
  33. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  34.      M       (input) INTEGER
  35.              The number of rows of the matrix A.  M >= 0.
  36.  
  37.      N       (input) INTEGER
  38.              The number of columns of the matrix A.  N >= 0.
  39.  
  40.      KL      (input) INTEGER
  41.              The number of subdiagonals within the band of A.  KL >= 0.
  42.  
  43.      KU      (input) INTEGER
  44.              The number of superdiagonals within the band of A.  KU >= 0.
  45.  
  46.      AB      (input/output) COMPLEX*16 array, dimension (LDAB,N)
  47.              On entry, the matrix A in band storage, in rows 1 to KL+KU+1.
  48.              The j-th column of A is stored in the j-th column of the array AB
  49.              as follows:  AB(ku+1+i-j,j) = A(i,j) for max(1,j-
  50.              ku)<=i<=min(m,j+kl)
  51.  
  52.              On exit, the equilibrated matrix, in the same storage format as
  53.              A.  See EQUED for the form of the equilibrated matrix.
  54.  
  55.      LDAB    (input) INTEGER
  56.              The leading dimension of the array AB.  LDA >= KL+KU+1.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZLLLLAAAAQQQQGGGGBBBB((((3333FFFF))))                                                          ZZZZLLLLAAAAQQQQGGGGBBBB((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      R       (output) DOUBLE PRECISION array, dimension (M)
  75.              The row scale factors for A.
  76.  
  77.      C       (output) DOUBLE PRECISION array, dimension (N)
  78.              The column scale factors for A.
  79.  
  80.      ROWCND  (output) DOUBLE PRECISION
  81.              Ratio of the smallest R(i) to the largest R(i).
  82.  
  83.      COLCND  (output) DOUBLE PRECISION
  84.              Ratio of the smallest C(i) to the largest C(i).
  85.  
  86.      AMAX    (input) DOUBLE PRECISION
  87.              Absolute value of largest matrix entry.
  88.  
  89.      EQUED   (output) CHARACTER*1
  90.              Specifies the form of equilibration that was done.  = 'N':  No
  91.              equilibration
  92.              = 'R':  Row equilibration, i.e., A has been premultiplied by
  93.              diag(R).  = 'C':  Column equilibration, i.e., A has been
  94.              postmultiplied by diag(C).  = 'B':  Both row and column
  95.              equilibration, i.e., A has been replaced by diag(R) * A *
  96.              diag(C).
  97.  
  98. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  99.      THRESH is a threshold value used to decide if row or column scaling
  100.      should be done based on the ratio of the row or column scaling factors.
  101.      If ROWCND < THRESH, row scaling is done, and if COLCND < THRESH, column
  102.      scaling is done.
  103.  
  104.      LARGE and SMALL are threshold values used to decide if row scaling should
  105.      be done based on the absolute size of the largest matrix element.  If
  106.      AMAX > LARGE or AMAX < SMALL, row scaling is done.
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.